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Abstract 

We reexamine the popular belief that a telluric planet or satellite on an eccentric orbit 
can, outside a spin-orbit resonance, be captured in a quasi-static tidal equilibrium called 
pseudosynchronous rotation. The existence of such configurations was deduced from 
oversimplified tidal models assuming either a constant tidal torque or a torque linear in the 
tidal frequency. A more accurate treatment requires that the torque be decomposed into 
the Darwin-Kaula series over the tidal modes, and that this decomposition be combined 
with a realistic choice of rheological properties of the mantle, which we choose to be a 
combination of the Andrade model at ordinary frequencies and the Maxwell model at low 
frequencies. This development demonstrates that there exist no stable equilibrium states 
for solid planets and moons, other than spin-orbit resonances. 



1 Motivation 

The ongoing quest for extraterrestrial life has placed exoplanets and their properties into 
the forefront of scientific investigation. The trend has provided additional momentum to 
the development of a broad variety of techniques and approaches employed in the planetary 
sciences. For example, the recent revival of interest in mechanics of bodily tides is partly due 
to the importance of the planetary spin for the prospects of finding habitable worlds near other 
stars. 

Well-known examples of dynamical equilibria achieved via tidal coupling include our Moon, 
which is in a 1:1 spin-orbit resonance to the Earth, and Mercury which makes exactly three 
sidereal rotations over every two orbital revolutions around the Sun. Similar behaviour is 
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expected of the growing number of known super-earths - especially if their composition happens 
to be similar to that of the terrestrial planets of the Solar System, i.e., if they have massive 
solid or partially molten mantles of rocky minerals. 

Unfortunately, some of the published far-reaching conclusions about specific exoplanets are 
based on incomplete or ad hoc models which should never be used for solid materials, including 
those with partial melt. Both these models, introduced by Goldreich (1966) mainly for the ease 
of analytical treatment, predict quasi-static pseudosynchronous rotation states, with the planet 
being trapped in a slowly changing equilibrium state at a faster-than-synchronous rotation rate 
and a vanishing orbit-averaged tidal torque. Except in specific (very narrow) frequency bands, 
these models are incompatible with the rheological properties of realistic mantles and crusts. 
Analysis based on actual rheologies demonstrates the impossibility of pseudosynchronous ro- 
tation for homogeneous terrestrial objects. Whether this prohibition extends to planets and 
moons with internal or surface oceans remains an open issue and needs further research. 

2 The constant angular lag model 

A consistent linear theory of bodily tides is based on Fourier decomposition of the tide, with 
subsequent inclusion of the response at each separate mode. The ensuing level of complex- 
ity has tempted many to circumvent it by developing simpler approaches. Serving as good 
illustrations and reflecting some qualitative aspects of the tidal interaction, such models are 
not necessarily applicable for quantitative purposes (Efroimsky & Lainey 2007) and should 
certainly be eschewed when fine features of near-resonant dynamics are explored. 

2.1 The essence of the method 

One, often-used, toy model prescribes: 

(a) to set both the Love number k% and geometric lag e g frequency-independent; 

(b) to insert their values into the popular short formula 

% = ^GMj 2 ^ k 2 sm2e g (la) 

for the polar component of the torque wherewith a tide-raising perturber of mass Mi 
acts on a tidally perturbed body of radius R located at a distance r , the obliquity i 
assumed small. 

(c) to combine formula (Hap with the assumption that the angle e g stays constant while 
the tide-raising perturber stays on one side of the bulge (in the sense of the directions 
as seen from the perturbed body's centre). 

For a nonzero eccentricity e, and in a sufficient proximity of the 1:1 resonance, the relative 
orientation of the perturber and the bulge changes twice over an orbital period. Hence, within 
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this model, the angle e g is set, by hand, to change its sign abruptly twice in a cycle, while 
keeping its magnitude | e g | fixed. Therefore the model can be written down as 

% = -GM? 1 ^ k 2 sm2\e g \ Sgn (z> - 9) , (lb) 

9 and 9 being the perturbed body's sidereal angle and spin rate, and v being the true anomaly. 

Historically, the method dates back to the paper by MacDonald (1964) and therefore is 
often referred to as the MacDonald torque (e.g., Touma & Wisdom 1994, section 2.7.1). 

The approach is also called the constant angular lag model or the constant tidal torque 
model, both names being somewhat misleading. Indeed, in the vicinity of the 1:1 resonance 
the sign of the lag (and the torque) is set positive or negative, when the bulge falls behind or 
advances relative to the direction towards the perturber. So both quantities change their sign 
twice over a period - a circumstance that makes the term constant inappropriate. Furthermore, 
the torque depends upon the distance and is always evolving in time unless the orbit is circular. 

The abrupt switch of the sign of the torque ([I]), with its magnitude staying unchanged, is 
quite a contrived assertion^ which by itself indicates that the model is unphysical. A deeper, 
mathematical, objection will be brought up in subsection 12.41 

Saying goodbye to the constant-angular-lag (constant-torque) model will not be easy, be- 
cause it has been a textbook standard for nearly half a century. Given the attractive simplicity 
of the model, one will always be tempted to enquire if perhaps it would still be producing at 
least qualitative results of some value. To see that it would not, we shall have to scrutinise the 
principal outcomes of the model. 

Specifically, we shall be interested in what happens when the tidal torque acts on a body 
of no permanent triaxiality. The perturber's orbit is set to lie in the equatorial plane of the 
perturbed body, i.e., the obliquity is set nil. Two special situations of interest emerge here. 
One is the case of exact synchronism, the other being the case of a vanishing average tidal 
torque. Both settings were explored by Goldreich (1966) whose results are explained in detail 
by Murray & Dermott (1999). 

2.2 The synchronous spin case 

Suppose the tidally-perturbed body on an elliptic (e / 0) orbit is caught into the 1:1 spin- 
orbit resonance: n = 9 . Then, as explained in subsection 5.2 of Murray & Dermott (1999), 
the angular motion rate 9 exceeds n over exactly one half of the orbital time period, and 
falls short of n during the other half of the period. Correspondingly, the tidal torque T z is 
positive (accelerating) through the former half of the period, and is negative (decelerating) 
through the latter half. The instantaneous tidal torque is proportional to a negative power of 
the instantaneous distance r between the bodies. As depicted in Figure 5.3 in Ibid, when the 
disturbed body's angular motion is faster than the mean motion, the bodies are closer, so the 
positive (spinning up) tidal torque is larger in absolute value than the negative torque for the 
other half of the period. Thus the resultant orbit-average torque ( T z ) is positive, and the net 
effect is to accelerate the tidally perturbed body's spin. (Recall that the undisturbed body is 
assumed to be spherical or oblate, so the tidal torque is the only one coming into play.) 

1 Stated alternatively, if we represent (a/r) 6 as a series of Fourier harmonics cos(jAA), j = 0, 1, . . ., where 
a is the scmimajor axis and A4 is the mean anomaly, we shall have to accept that the cos(A / I) tidal mode 
generates a positive (accelerating) torque for Ai E [— tt/2, +tt/2], abruptly switching to a negative value for 
M G [7r/2,37r/2] . 
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2.3 The case of vanishing tidal torque 



The second important application of the constant angular lag model is the situation where 
the orbit-average tidal torque vanishes: (71) = . Once again, the body being spherical or 
oblate, no triaxiality-caused torque shows up. So vanishing of the average tidal torque entails 
a dynamical equilibrium: the tidally disturbed body keeps rotating at a steady spin rate. A 
calculation of this rate, borrowed from Goldreich (1966), is presented in Murray & Dermott 
(1999) and is often cited in the literature. According to that development, the equilibrium is 
achieved, for a zero obliquity, at the spin rate of 

^="(1 + y^j , (2) 

e being the eccentricity. At first glance, the result looks unassailable. Indeed, for 6 = n , the 
bulge is lagging behind the central line over one half of the time period (around the periastron), 
so the torque accelerates the rotation. Over the other half of the period, the torque decelerates, 
being weaker due to a larger distance. So the state 9 = n looks unstable, as the overall average 
torque seems to be accelerating. 

It is however well known that our own Moon is not staying in this pseudosynchronous 
regime (which would be 3% faster than the synchronous rotation wherein the Moon is presently 
locked). Murray & Dermott (1999) point at the lunar quadrupole moment as the reason why 
the Moon is not pseudosynchronous. A deeper reason though lies in the constant geometric 
lag model being genuinely flawed, and in the entire calculation leading to (j2J) being invalid. 

2.4 A major objection against the constant-angular-lag model 

As well known, the generic expression for the tidal amendment to the perturbed body's poten- 
tial is furnished by a Fourier series developed by Kaula (1964). We term it the Darwin-Kaula 
expansion, as a partial sum of that series was written down by Darwin (1879). Accordingly, 
the generic expression of the tidal torque also must look as an infinite series. The series re- 
mains infinite even if we include into it only the degree-2 terms, i.e., those proportional to the 
quadrupole Love number hi . The very fact that the expansion for the torque can be wrapped 
into a short and neat form ([T]) is an indicator of some extra, very special assumption being 
involved. 

As was pointed out in Williams & Efroimsky (2012), such an assumption indeed is present in 
the constant-angular-lag (constant-torque) model, though this assumption is never stipulated 
explicitly. The situation is elucidated in all detail in the paper by Efroimsky & Makarov (2012) 
to which we refer the reader. Here we shall provide only a brief summary. 

As explained in Williams & Efroimsky (2012), the concise expression ((T]) is equivalent to 
the full Darwin-Kaula expansion for the potential, only if the following assumptions are made: 

• In all terms of the Darwin-Kaula series for the tidal amendment to the potential of the 
perturbed body, i.e, for all Fourier tidal modes ui mpq , the time lags are endowed with 
the same, frequency-independent value At . 

• The obliquity is set small. 
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• Only the terms with I = m = 2 are retained, and the Love number hi entering these 
terms is assumed frequency- independent. Here the degree / and the order m are the 
first two integers of the four-number set Impq used to number the Fourier modes showing 
up in the Darwin-Kaula expansion. 

The so-processed Darwin-Kaula series for the tidal potential becomes equivalent to a concise 
expression (eqn. 16b in Ibid.) wherefrom our expression ([1]) for the torque ensues. 

Alternatively, the above three assumptions could be applied directly to the Darwin-Kaula 
expansion for the tidal torque. Once again, the outcome would be the above expression ([I]) for 
the torque. This is demonstrated in Efroimsky & Makarov (2012, eqn. 29). 

Under the three assumptions, the instantaneous geometric lag angle turns out to be @ 

e g = (z> - 9) At , (3) 

At being the frequency-independent time lag, v being the true anomaly of the perturber, 
and 9 being the sidereal angle of the tidally perturbed body. From this expression, it is 
straightforward that the validity of formula ([T]) is incompatible with the geometric lag being 
constant. Indeed, the road to is paved with the aforementioned three assumptions, one of 
which being that of a constant At. As can be observed from ([3]), the latter is incompatible 
with the geometric lag being constant, unless the eccentricity is nil. H 

On all these grounds, the constant-geometric-lag (constant-torque) model should be dis- 
carded as such. 

3 Pseudosynchronism in the constant time lag model 

As distinct from the constant geometric lag approach, the constant time lag model sets the 
time delay At independent of the tidal mode frequency. Pioneered by Darwin (1879), this 
assumption was a part of numerous works, e.g., Hut (1981), Eggleton et al. (1998). The 
assumption was also the base for one of the two models considered by Goldreich & Peale 
(1966, eqn. 23), the other model addressed in that paper being the afore-discarded constant 
geometric lag method. 

The constant time lag model is unique, in that it makes the full Darwin-Kaula expansion for 
the tidal potential (or torque) equivalent to a much shorter and simpler expression. In regard 
to the tidal potential, this is the equivalence of the full series (16) and a simpler expression 
(15) in our preceding paper Efroimsky & Makarov (2012). In application to the torque, this is 
the equivalence of the appropriate full series to the simpler expression (23) in Ibid. 

2 The geometric angular lag e g is not to be confused with the instantaneous phase lag (or longitudinal lag) 

e ph = 2(u-0)At = 2e g 

sometimes used in the literature (Efroimsky & Williams 2009, Williams & Efroimsky 2012). 

3 As a last resort, one can suggest (a) to tune the time dependence of At so that the lag angle e g in ^ 
stays constant in time, and (b) to assume that the time lags at all Fourier tidal modes are equal to the so 
specially evolving At . This would imply that the lagging properties of the material at all frequencies are being 
tuned in a fine manner, continuously and simultaneously, so that At adjusts its evolution rate, to stay inverse 
to v — at any instant of time. The rate of change of v — 6 being defined by the orbital parameters, existence 
of such a rhcology in nature looks impossible. 
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If we agree to limit our approximation to the lowest degree and order, I = m = 2 , the 
aforementioned simpler expressions become especially neat and concise. They read as 
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U{r) = - ^GM 1 k 2 — cos( 2 [y - 9) At ) = ^ G M 1 k 2 ^ cos( 2 | v - 9 | At ) , (4) 
for the potential, and as 

T z = - G Ml k 2 sin( 2 (is - 9) At) = - G k 2 ^ sin( 2 | v - 9 | At ) Sgn (z> - 9) , (5) 

for the polar torque. A detailed derivation of fl4]) can be found in Williams & Efroimsky (2012), 
while derivation of (J5J) is offered in Efroimsky & Makarov (2012). 

To average the torque, it is instrumental to insert into flS]) the distance r expressed through 
the semimajor axis a , eccentricity e , and true anomaly v , and to integrate over the orbital 
cycle. The procedure gets simplified greatly for small lags, when one can substitute the sine 
with its argument. Then the calculation (presented in detail in the Appendix to Williams &; 
Efroimsky 2012) renders: 



X) oc 



f e 2 +f e 4 +|e 6 Q 1 + 3 e 2 + | e 4 



e 2f n (1 _ ^ 
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This expression was obtained by Eggleton et al. (1998), though its equivalent was present in 
an earlier paper by Hut (1981, eqn. 11). In a somewhat disguised form, this expression can be 
found in a much earlier work by Goldreich & Peale (1966, eqn. 24). 

We find readily that the equilibrium (zero average tidal torque) is achieved at 
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Note that the pseudosynchronous rate of rotation depends only on the mean motion and 
eccentricity. This enables us to solve equation (J7J) with respect to e , for a fixed dimensionless 
spin rate 9/n . The outcome will be the equilibrium eccentricity e equ , i.e., the eccentricity that 
ensures the vanishing of the average torque at a certain value of 9/n. In Figures [1] and El 
the equilibrium eccentricity e equ is depicted as a function of 9/n. For the model leading to 
expression ([6]) for the torque, this is a monotonically rising curve. 

The curve divides the plane into two parts corresponding to the two opposite signs of the 
average polar torque ( % ) . While ( % ) is positive (accelerating) everywhere above the curve, 
is stays negative (decelerating) everywhere below the curve. Indeed, if we fix the eccentricity 
and make 9/n very large, this will guarantee us that we get into the lower right part of 
the picture, i.e., below the rising curve. In this situation, i.e., for a fixed eccentricity and a 
sufficiently swift spin, the second term of the torque ([6]) must be leading, wherefore the torque 
must be negative, i.e., despinning. Similarly, by fixing the eccentricity and making the spin 
rate very small, we ensure getting into the upper left part of the picture, and also ensure that 
the first term in ([6]) is leading, so the torque is positive, i.e., accelerating the spin. Since the 
smoothly rising curve0 e equ (9/n) corresponds to a zero (%} , it is impossible to change the 
sign of ( % ) without crossing the curve. 
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1.0 1.2 1.4 1.6 1.8 2.0 



rate of rotation 0/n 

Figure 1: The equilibrium eccentricity (one corresponding to a vanishing average tidal 
torque) depicted against the dimensionless spin rate 0/n . Calculations were made for 
a tidally perturbed rotating body with parameters of the Moon, as shown in Table [U 

The monotonously rising curve corresponds to the linear torque (constant time 
lag) model. The jagged dotted line corresponds to a realistic rheology intro- 

duced in Section 01 Both functions were computed with a step of 0.01 in 0/n. 

In both cases, the resulting curve divides the plane into two parts corresponding to the 
two opposite signs of the average polar torque ( T z ) ■ While (T z ) is positive (accelerat- 
ing) everywhere above the curve, it stays negative (decelerating) everywhere below the curve. 

The small arrows indicate the action of the tidal torque upon small perturbations of the spin rate 
away from an equilibrium state. For the constant time lag model, the torque is restoring, and the 
equilibrium is stable. In the case of realistic rheology, though, the emerging nonzero torque drives 
the rotator away from the stable spin. 



7 



T 




1.0 1.2 1.4 1.6 1.8 2.0 

rate of rotation 6/n 

Figure 2: Equilibrium eccentricities of a zero secular tidal torque acting on a tidally perturbed 
super-Earth, depicted against the dimensionless spin rate 6/n. Parameters of the super-Earth 
are given in Table [IJ and are consistent with those chosen for GJ581d in Makarov et al. (2012). 

The monotonously rising curve represents the prediction of the constant time lag model. 

The jigsaw dotted curve illustrates the prediction of the realistic rheological model described in 
Section HI The function was computed for a grid of points at a step of 0.01 in 6/n . 
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Within the constant time lag model, the function e equ (8/n) being a smoothly rising curve 
explains the emergence of pseudosynchronism. To see this, consider a point on this curve, 
corresponding to a certain pseudosynchronous state. A small perturbation in 9/n makes the 
tidally perturbed body rotate either faster or slower than the pseudosyncronous rate, and a 
nonzero tidal torque emerges. Illustrated by the two counter directed short arrows on the plot 
in Figured] the tidal torque is restoring, in that its action is opposite to the sign of perturbation. 
Thus, the tidal torque will return the perturbed rotator to the equilibrium state, i.e., to the 
initial position on the curve. So the equilibrium is stable. % 

The constant time lag model ignores the important contribution of rigidity (Segatz et al. 
1988) and inelasticity (Karato and Spetzler 1990) into the tidal response of Earth-like planets. 
As a result, the model is incapable to account correctly for creep. As will be discussed in 
the following section, the above derivation of quasi-stable pseudosynchronism, from the linear 
torque model, is inapplicable to Earth-like planets with rigid mantles. However, in the viscous 
limit, this model may still be applicable to celestial bodies that do not have appreciable rigidity 
or inelasticity, such as gaseous planets and stars. Observations of binary stars, especially of 
short-period active stars on eccentric orbit, hold the best prospect of proving or disproving the 
linear torque model for this type of objects (Ferraz-Mello 2012, Torres et al. 2010). 

The spin rate of active stars can be inferred from the characteristic periods of photometric 
variations caused by the passage of large spots or groups of spots across the visible disk of the 
star. The orbital period and the eccentricity are determined from spectroscopic radial velocity 
measurements. We find somewhat conflicting evidence for the existence of pseudosychronism 
in binary stars. Some stars with considerable eccentricities appear to have pseudosynchronous 
rotation (Hall 1986, Fekel et al. 1998), which is consistent with the original prediction by 
Hut (1981). Other stars clearly rotate faster or slower than the predicted rate (Fekel et al. 
1993, Strassmeier et al. 2011). Even more puzzling, a significant number of tight binary 
systems have been found on circularised orbits, albeit spinning clearly asynchronously. This 
fact comes into contradiction with one of the important predictions of the linear torque theory, 
the one that synchronisation (or pseudosyncronisation) of rotation is achieved much sooner 
than circularisation (Eggleton 2011). Thus the impression created by the current body of 
observations is that the constant time lag model is, at least, not universally applicable to stars. 
This should not come as a surprise, because there exist theoretical indications that stars may 
have magnetic rigidity (Williams 2004, 2005, 2006; Ogilvie 2008, Garaud et al. 2010). E 

4 Here and hereafter, the symbol e equ (9/n) implies e equ as a function of the ratio 9/n. This is not a 
product of e equ and 9/n. 

5 Being stable, the equilibrium is gwasi-static, in the following sense. As the tidal dissipation goes on, the 
process of despinning continues. The argument 9/n slowly decreases, and so does the appropriate value of the 
equilibrium eccentricity. After a perturbation in e from an equilibrium state gives birth to a torque, the torque 
corrects swiftly the spin rate in such a way that the rotator returns to an equilibrium state. However, the 
equilibrium state itself is evolving slowly. In the case of a two-body problem, this evolution is always directed 
towards the configuration with 9/n — 1 and e = . For a viscous body, this was proven by Hut (1981). For 
a broader class of viscoelastic rheologies, the proof was offered by Bambusi and Haus (2012). 

6 Another deviation from the purely viscous (constant time lag) model can be caused by the so-called 
A— effect responsible for differential rotation (Kapyla & Brandenburg 2008, Kichatinov 2005, Riidiger 1989). 
Turbulent convection generates an extra stress called Reynolds stress. While in a non-rotating convection zone 
this stress can be described as an addition to the viscosity, this can no longer be done when the rotation period 
becomes comparable to or shorter than the convective turnover time. In that situation, a non-viscous input, 
the so-called A-effect, shows up. In its presence, the stress tensor in the stellar material will no longer be 
proportional to the time derivative of the strain tensor, but will contain terms proportional directly to velocity. 
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Table 1: Parameters of the tidal model. 



Name Description Units Values 

Moon super-Earth 

(GJ581d) 

£ moment of inertia coefficient 2/5 2/5 

R radius of planet m 1.737 ■ 10 6 1.083 • 10 7 

M 2 mass of the perturbed body kg 7.3477 ■ 10 22 4.23 ■ 10 25 

Mi mass of the perturbing body kg 5.97 • 10 24 6.17 ■ 10 29 

a semimajor axis m 3.84399 • 10 8 3.3 ■ 10 10 

n .. mean motion, i.e. 27r/P or b yr -1 84 34.25 

e orbital eccentricity 0.0549 0.27 

(B-A)/C triaxiality 2.278 • 10" 4 5 • 10" 5 

G gravitational constant m 3 kg" 1 yr" 2 66468 66468 

T^j Maxwell time yr 5 50 

fi . unrelaxed rigidity modulus Pa 0.8 ■ 10 11 0.8 • 10 11 

a the Andrade parameter 0.2 0.2 



4 Equilibrium torques for Earth-like planets and moons 

To build a consistent theory of bodily tides, one has, first, to decompose the tide into a 
Fourier series and, second, to attribute to each Fourier component its own phase delay and 
magnitude decrease (the latter being expressed by the Love number appropriate to the said 
Fourier mode). Development of the decomposition technique was started by Darwin (1879) 
and accomplished in full by Kaula (1964). Attribution of phase delays and Love number values 
to the Fourier modes took much longer time, because of the necessity to explore rheological 
properties of the mantle at various frequencies. This exploration, by both seismological and 
geodetic methods, has been going on intensively through the past dozens of years. Merger 
of the Darwin-Kaula decomposition technique with the results from solid-Earth rheology is 
explained in Efroimsky (2012 a). The paper relied on a combined rheological model (Andrade 
at higher frequencies, Maxwell at lower frequencies), because of this model's ability to best 
match laboratory experiments and both seismic and geodetic measurements of dissipation over 
a range of frequencies in the solid Earth0 The merger of the Darwin-Kaula expansion with the 
combined rheological model, worked out in Ibid., has been used to predict spin-orbit resonances 
of a Mercury analogue having a constant eccentricity and a zero obliquity (Makarov 2012) and 
tidal properties of super-Earths (Efroimsky 2012 b). 

Thus the purely viscous model falls apart, and a frequency- independent time lag is no longer an option. 

7 Motivation for the combined model stems from the mantle being predominantly viscoelastic at frequencies 
below some threshold, and predominantly inelastic at frequencies above it. As explained in Karato & Spetzler 
(1990), dissipation above the threshold is dominated by defect unpinning (see also Miguel et al. 2002). When the 
frequency descends below the threshold, the effectiveness of this mechanism declines, because the Andrade term 
in the expression for the complex compliance decreases exponentially. The response of the mantle approaches 
that of the Maxwell body. So slow processes (like the postglacial rebound) are viscoelastic. 

For Earth's mantle, the threshold frequency is of the order of 1 yr" 1 . Its value, though, is exponentially 
sensitive to the temperature and therefore may be very different for exoplanets and exomoons. 
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4.1 Expression for the tidal torque 

It is explained in Appendix A that the average polar component of the tidal torque can be 
approximated with the following expression, provided that (1) the obliquity is small, and (2) 
the perturbed body and the perturber are not too close to one another (so only the terms with 
degree-2 Love number are important) : 



^ ' q=-X 

l^T (f) 5 ^ G l q ^ k ^ m ) sin I^K 2 og)l Sgn(c 220g ) + O(e 8 e ) + O(z 2 e ) . (81 

^ ' q=— 1 

This is the polar (orthogonal to the equator) component of the torque wherewith the tidally- 
perturbed body is acted upon by the perturber. The angular brackets denote orbital averaging, 
G stands for the Newton gravitational constant, Mi signifies the mass of the perturber (the 
star, if the perturbed body is its planet; or the planet, if the perturbed body is a satellite), 
a is the semimajor axis, while R is the radius of the tidally perturbed body. The degree-2 
dynamical Love number fc 2 and the phase lag e 2 are functions of the Fourier tidal mode 

^ 22 og = ( 2 + 9) n ~ 2 * " ( 9 ) 

While the dynamical Love numbers k 2 (u 220 q) are positive definite, the sign of each phase lag 
e 2 (w 220O ) coincides with that of the Fourier mode be understood from formula 

fl20l in Appendix A. It is for this reason that the products k 2 (u 220 g) sin e 2 (w 220 g) emerging in 
expression fl8a|) are rewritten in expression fl8bl) as k 2 (u 220 q) sin | e 2 (w 220 g) | Sgn (w 220 ^) ■ 

A generic expression for the torque implies summation over the four integer indices Impq 
serving to number the Fourier tidal modes u) lmpq entering the spectrum - see Appendix A 
below. The terms of that series depend also upon the obliquity. As the Impq term contains 
a factor (R/a) 2l+1 , it is often sufficient to keep only the degree-2 terms (1 = 2). In this case, 
with an extra assumption of small obliquity, it is enough to limit the summation to the terms 
with m = 2 , p = . This renders expression ([8]). 

While the full expression for the torque implies summation over all integer values of q , 
numerical tests demonstrate that, for eccentricities not exceeding ~ 0.3 , it is enough to 
take into account the terms up to e 7 , inclusive. This would require summation from q = 
— 7 through q = 7 . However, the values of the numerical factors entering the eccentricity 
polynomials C 20 g( e ) are such that in practice it turns out to be sufficient to include only the 
terms with q varying from — 1 through 7 . 

4.2 The tidal torque and the equilibrium eccentricity as functions 
of the spin rate 

Expression (Q makes each product /c 2 sin e 2 a function of the planetary spin rate 9 : 
^Kaog) sin | e 2 (u 220q ) | Sgn ( u 220q ) 

= k 2 {2(n - 6) + qn) sin | e 2 ( 2{n - 6) + q n ) \ Sgn ( 2{n - 6) + q n ) . (10) 
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Consequently, the entire sum (JSJ) can be treated as a function of 9 . The mean motion and 
eccentricity will play the role of parameters whose evolution is much slower than that of 9 . 
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Figure 3: Angular acceleration due to the secular tidal torque ([8]), in the vicinity of the 2:1 resonance. 
The dotted kink is the q = 2 term which is an odd function when centered at 9/n = 1 + q/2 = 2. 
The solid line renders the total torque ([8]), i.e., a sum of the q = 2 kink and the bias comprised by 
the terms with q ^ 2. Near the resonance, the bias is a slowly changing function of 9/n and can be 
approximated with a constant. The q = 2 kink resides on the right slope of a more powerful q = 1 
kink which is centered at 9/n = 1 and dominates the bias. So the q = 2 kink is shifted downwards 
and goes through nil a tiny bit to the left of the resonance. 

The figure is borrowed from our work Makarov et al. (2012) devoted to the superearth GJ581d. 



Each product k 2 sin e 2 is an odd function of the tidal mode w 22 og and has the shape of 
a kink centered around oJrmq = . When we employ relation ( TTUT) to write these products 
as functions of the spin rate, the new functions will still be kinks, though centered around 
9 = n (1 + q/2) . In Figure [31 the dotted line depicts the product 
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The kink shape of the fc 2 sin e 2 products is determined by the rheological properties of the 
planet and its self-gravitation (see Appendix B for details and references). A kink transcends 
nil and changes sign continuously as the spin rate goes through the appropriate resonance. 
In the sum ([8]), the kink-shaped products stand with multipliers ^^^(e) . So the overall 

tidal torque (IE]), as a function of 9 , is a superposition of many kinks having different magnitudes 
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and centered at different resonances (nine kinks, if we sum over q = — 1, ... 7). The ensuing 

curve will cross the horizontal axis in points extremely close to the resonances 9 = n (1 + q/2) , 
but not exactly in these resonances — like the solid line in Figure El 

4.3 The physical meaning of the kink 

The physical forcing frequencies in the mantle, %22og > are absolute values of the Fourier modes: 

^220g = I U 220q I • (12) 

The dynamical Love number is an even function of the tidal mode ui 2 2oq > wnne the phase lag 
is an odd function. For this reason, the product k 2 (u 220 q) sin £i{u) 220 q) can ^ e rewritten as a 
function of the physical frequency x , multiplied by the sign of the appropriate Fourier mode: 

M^jaog) sme zK 2 og) = k 2(X 2 2oq) sin l e 2(x 22 og)l S § n K 2 og) ■ ( 13 ) 

Outside the inter-peak interval (i.e., at frequencies that are not too low), the positive definite 
quantity k 2 (x) sin|e2(x)| is decreasing monotonously with increase of the frequency \- 
This happens for two reasons. One, intuitively obvious, reason is that the dynamical Love 
number decreases at higher frequencies, because materials are inertial, and it is getting harder 
for their shape to keep up with the varying stress as the frequency goes up. Less obvious is the 
circumstance that the sine of phase lag (i.e., the inverse tidal quality factors), too, decreases as 
the frequency growsj^l Supported by a mighty volume of seismological, geodetic, and laboratory 
data, this behaviour may look counterintuitive because this is not what one would expect from a 
viscous fluid. The fact however is that at physically interesting frequencies the mantle behaves 
itself not as a viscous or a Kelvin- Vogt body but as an Andrade body dissipation wherein 
obeys the law sine oc x~ a , with a ~ 0.14 — 0.4 for most solids (Efroimsky 2012a,b). 

Finally, the steep (but still continuous) near-resonant jumps connecting the peaks of the 
kink are explained by the fact that at those locations self-gravitation "beats" rheology (Ibid.). 

The kink shape of k 2 sin e 2 (generally, of k\ sin ei ) entails somewhat counterintuitive con- 
sequences for the phase lag and the geometric lag angle. Consider the principal, semidiurnal 
bulge. Its phase lag and the geometric lag angle are 

e 22oo = ^2200 ^2200 = 2 (n — 6>) At 2200 (14) 

and 

^2200 = 2 I £ 2200 I = \ U ~ Q\ ^2200 ' (^) 

8 This quantity is often denoted as k-2/Q , though notation &2 Qi would be more appropriate. Tidal quality 
factors are not identical to the seismic quality factor, the difference becoming crucial at low frequencies. 

9 To illustrate the decrease of the Love number, imagine that we dip a spoon into a bowl of honey, and apply 
to the spoon an oscillating force of a fixed amplitude. Naturally, the amplitude of motion of the spoon will be 
larger for lower frequencies. Sadly, this simple example will not help us to illustrate how the phase decreases 
with the growth of frequency. Naively, one might expect an anti-phase response at high frequencies, like in the 
case of a damped driven harmonic oscillator. This regime would indeed be taking place, had the mantle obeyed 
a constant time lag law. However, real minerals behave differently, so our parallels with a viscously damped 
oscillator or a viscous liquid have their limitations. 
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Were a planet composed of a material with the time lag At 2200 insensitive to the value of 
the principal forcing frequency X2200 = ^ \ n ~ ® \ ■> ^ ne geometric lag angle would be larger 
for a higher value of this frequency. This indeed is what one would, intuitively, expect: the 
higher the difference between 9 and n the larger the angle. However, for a realistic rheology, 
an increase of the <5 2 2oo angle due to an increase in 2\n — 9\ will take place only within an 
extremely close proximity of the the 1 : 1 resonance. Stepping beyond the kink's peak, we 
shall find that an increase in 2\n — 6\ will be accompanied by such a decrease in the time 
lag At 2200 that the product of these two quantities will, overall, be decreasing with growing 
frequency. So both the phase lag and the geometric lag angle will become smaller. 

4.4 Instability of pseudosynchronous rotation. Physical interpreta- 
tion 

Since the mode-dependence of the products (TTUT) follows from the rheological properties of the 
mantle and from self-gravitation of the planet, these products turn out to be functions not 
only of the tidal mode u 2 2oq > but a ^ so °f ^he parameters defining the size and rheology of the 
perturbed body. These parameters (presented in Table 1) are the body's radius R and mass 
M 2 , as well as its unrelaxed rigidity fi , Maxwell time tm , inelastic (Andrade) time T4 , and 
the Andrade parameter a . 

For a given selection of these parameters' values and a fixed spin rate 9 , there is a unique 
eccentricity e equ at which ( T z ) = . The dependence of e equ on the relative rate of rotation 

9/n can be found for a grid of points by numerically determining the roots of (T z ) in 

e . To accomplish this, we recall that each term in equation (JSJ) is a polynomial in e and, 
therefore, so is the entire right-hand side of (JSJ). Computational search of the roots was carried 
out for two sets of parameters listed in Table 1, one representing the Moon orbiting the Earth, 
and the other a hypothetical super-Earth orbiting a solar analog. The choice of parameters is 
intended to represent the range of applicability of the model. The resulting dependencies of 
e C qu upon 9/n are presented in Figures [1] and |2j 

The jigsaw shape of the found dependencies e cqn (9/n) is remarkably different from the 
predictions of the linear torque model. The curves for the Moon (Fig. [TJ) and the super-Earth 
(Fig. [2]) are monotonously descending with a rising rate of rotation, everywhere between the 
spin-orbit resonances. This has profound implications for the character of equilibrium at a 
zero tidal torque. Consider an arbitrary point on the downhill portion of the curve, e.g., the 
one in Fig. [H from which two short arrows of opposite direction are drawn. A perturbation 
in 9 away from this point, whether spinning the planet up or slowing it down, will cause a 
nonzero tidal torque acting in the same direction - as indicated by the direction of the arrows. 
Thereby, the tidal equilibrium achieved at e equ between the resonances is inherently unstable. 
Similarly, a perturbation in e will make the planet diverge from the curve of zero torque rather 
than return to it. Hence, the states of zero tidal torque at non-resonant spin rates are transient 
by nature. 

In a narrow vicinity of each spin-orbit resonance, e equ takes a rapid upward increase. 
Computation of the roots of Equation [S] is numerically difficult in these areas because of a very 
large gradient of the curve. The slope of the segments at the resonances is positive; therefore, 
there exists a stable equilibrium just as in the case of the linear torque model. A deviation of 
the spin rate from a resonant value gives rise to a nonzero restoring torque, as indicated by a 
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pair of inward arrows in Figured] at the resonance 3:2. 

To obtain a physical explanation of the unstable nature of the pseudosynchronous rotation, 
recall two circumstances. First, each term of the sum (jSj) has the shape of a kink. Second, 
each such kink, as a function of 9 , is increasing monotonously everywhere except near the 
appropriate resonance, as in Figure El For this reason, an infinitesimal increase in 9/n furnishes 
an infinitesimal increase in the tidal torque (not necessarily in its absolute value) . This happens 
for an arbitrary value of e and for the values of 9/n outside the narrow resonances. Specifically, 
for e = e equ , the torque is zero and acquires a positive value, which leads to further spin-up. 
The spin- up continues until 9/n stumbles into a resonance. (Resonances are depicted with 
near-vertical segments of the dotted curves in Figures [T] and [2J) 

4.5 On the choice of the values for physical parameters 

Figures [1] and [2] reveal, in comparison, that the structure of the equilibrium tidal torque is 
similar for a wide range of values of planetary parameters. The values employed to build these 
plots are shown in Table [TJ The rapid jumps of e equ at resonances and the smooth descents 
between the resonances are characteristic of small moons and large planets likewise. The values 
of e equ at resonances appear to be the same for the Moon and the model super-Earth. The 
most noticeable difference is in the amplitude of the resonance jumps, which is significantly 
higher for the Moon. By experimenting with the input parameters, we found out that this 
amplitude is sensitive mainly to the Maxwell time tm , which differs by an order of magnitude 
between our model bodies. 

The choice of such a small Maxwell time for the Moon, only 5 years, is justified by the 
likely presence of a high percentage of partial melt in the lower lunar mantle. The presence of 
partial melt follows from the modeling carried out by Weber et al. (2011) and also from an 
earlier study by Nakamura et al. (1974). There exists data pointing at the possibility of the 
lunar Maxwell time being of the order of months. 

The choice of the parameters of the super-Earth was consistent with that made in Makarov 
et al. (2012) for the planet GJ581d. 

The Arrhenius law requires that planets and moons with hotter interiors have lower viscosity 
of mantles and, thus, have shorter Maxwell times. So we surmise that the secular tidal torque 
for such objects should be relatively more efficient in capturing at higher spin-orbit resonances. 

5 Resonant rotation of axisymmetric bodies 

We have determined that a stable spin-orbit equilibrium is achieved at spin rates where the 
value of the secular polar tidal torque is zero, while the derivative of the equilibrium eccentricity 
with respect to the spin rate, de equ / d{9 /n) , is negative. With the realistic tidal model 

1 The smallness of the lunar tm ensues from the unexpected frequency dependence of the lunar tidal 
Q factor. According to Williams et al. (2008), the tidal Q increases from ^29 at a month to ~ 35 at 
one year, a slope incompatible with the seismic properties of rocks which are expected to have a seismic Q 
decreasing with increase of the period. As explained in Efroimsky (2012a), this unexpected slope may have 
emerged due to the difference between the tidal and seismic Q at low frequencies. It is possible that the 
frequency range in which the lunar tides were studied could be close to or slightly left of the peak of the 
function ^2(^2200) sin | £2(0-12200) I • I n this case, Tm of the Moon should be about 0.1 year. Fortunately, the 
choice of the value of 7^/ does not influence considerably the jagged shape of the dependency e cqu (9/n) . The 
plot for tm = 0.1 yr turns out to be similar to the jagged plot for Tm = 0.1 yr depicted in Figure [TJ 
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discussed in Appendix B, this may happen only in the vicinity of spin-orbit resonances because 
the derivative of the torque is positive elsewhere. The secular torque (|HJ) and the angular 
acceleration caused by it, has "kinks" in the vicinity of spin-orbit commensurabilities 9/n = 
1 + q/2 , with an integer q of either sign. An example thereof is shown in Figure H] for 
the super-Earth model, within a segment of the spin rate around 9/n = 5/2. The kink is 
comprised by a local maximum below the resonance and a minimum above the resonance. 

To understand the plot in Figure HI recall that in the vicinity of each resonance q' , i.e., for 
9/n being close to 1 + q'/2 , the right-hand side of flS} can be decomposed into two parts. One 
part is the q = q' term. It is a kink-looking odd function of the tidal mode Co>220g , an d it goes 
through nil at exactly u)220q — . Due to OH]), this term can also be interpreted as a function 
of the spin rate, antisymmetric with respect to the resonance point 9/n = 1 + q'/2 where 
this term goes through nil. The second part, called bias, is the rest of the sum, i.e., the input 
of all the q ^ q' Fourier modes into the values assumed by the torque in the vicinity of the 
q = q' resonance. The bias can be negative or positive in value, depending on the eccentricity. 
For not too large eccentricities, it is usually negative. Being a very slowly changing function 
within the resonance interval, it can, to a good approximation, be assumed constant there. 

Having summed up all the terms in ([8]), and exploring the behaviour of this sum near 
9/n = 1 + q'/2, we see that the resulting curve does not cross the horizontal axis at 9/n = 
1 + q'/2. One can say that the bias slightly displaces the location of zeros. The zeros are 
located close to resonances but not exactly in resonances. 

In Figure Hh, the values of the overall torque (in fact, of the total angular acceleration 
proportional to the torque) are defined mostly by the q = 3 term which, in this vicinity, 
looks like an antisymmetric kink. However, the curve is shifted down due to the bias which is 
defined mainly by the right slope of the q = 1 kink located to the left. Since the right slope 
of the q = 1 kink is negative, the q = 3 kink is shifted down. As a result, the maximum 
secular torque barely rises above zero, and the curve has two zeros located to the left of the 
point 9/n = 5/2, close to the maximum of the kink. The interval of resonance is between 
9/n = 2.4998 and 2.5002. 

Figure Hb is a blow-up of Figure showing in greater detail the area around the maximum 
of the kink. The root of the function within the resonance interval is at 2.49985 rather than 
exactly 2.5 . For the same reason, the torque value is negative at 9/n = 2.5 . 

In the framework of this model, it is reasonable to define capture in resonance as an equi- 
librium state in which the body's average rotation rate stays between the values corresponding 
to these two peaks. This definition is adequate because, as we saw above, stable equilibrium 
is possible only on negative slopes of the angular acceleration (or tidal torque) as a function 
of the spin rate. 

When a triaxial planet is captured in a 5:2 resonance, its time-averaged spin rate is exactly 
2.5 n. Similarly, the Moon's spin rate, captured in synchronous rotation, is exactly In. 
However, we know that the time-averaged tidal torque is nonzero at this spin rate. Why does 
not the Moon accelerate? For a triaxial body, the nonzero secular tidal torque is compensated 
by a counteracting triaxiality-caused torque, through a small tilt of the average inertia axis 
with respect to the line connecting the centres of the bodies. A nonzero time-average tilt 
generates a secular triaxial torque. This mechanism of torque compensation obviously does 
not work for the rather hypothetical case of axisymmetric body, which would be subject to 
only tidal forces. Would the Moon be facing the Earth always with the same side if it were 
completely axisymmetric? First, we have to find out if capture in spin-orbit resonance is at all 
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Figure 4: Angular acceleration of the super-earth (Table 1) caused by the secular tidal torque in the 
vicinity of 5:2 resonance, (a) showing the entire resonance interval, and (b) showing in more detail 
the same curve in the area of the local maximum. 
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Figure 5: Capture of an axisymmetric super-earth (Table 1, but with (B—A)/C = 0) in 5:2 resonance. 
Note that the ultimate equilibrium spin rate is slightly less than 2.5 n. 



possible. The answer is yes, as long as the secular torque changes sign in the vicinity of that 
resonance 0. Even though the maximum torque in Figure H] barely rises above zero at spin 
rates below the resonance, it turns out to be sufficient for capture in 5:2. Figure \5\ displays the 
results of a numerical integration of the evolution of spin rate for the super-Earth model (Table 
1) at tm = 50 yr with an initial spin rate of 9(0) = 2.51 n . The planet is captured in about 
8500 yr, but the equilibrium spin rate at 2.49985 n is clearly below the resonance value. This 
value is consistent with the root of the secular torque within the resonance interval, Figure Hb. 
Thus, the equilibrium resonance state of an axisymmetric body is achieved at the spin rate 
where the secular tidal torque equals zero, as expected. 

11 Triaxial bodies can be captured in spin-orbit resonances even if the secular tidal torque is negative every- 
where in the vicinity of that resonance. 
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6 Word of caution 



As was demonstrated above, stability of pseudosynchronous spin hinges upon rheology. Being 
stable for the constant time lag model, the regime is expected to be transient for realistic 
mantles, because their rheology is expected to have one pronounced peak (not to count the 
opposite one at the negative value of the tidal mode) - see Figure [31 The situation will have 
to be re-analysed for bodies of complex structure (with surface or internal oceans), as well as 
for bodies of yet unexplored rheologies. Specifically, if it happens that somewhere in universe 
there exist bodies with not one but two pronounced peaks of k 2 sin e 2 at positive frequencies, 
the "ditch" between these peaks may, in principle, lead to emergence of a pseudosynchronous 
rotation state. Such a peak may emerge at the boundary of two frequency bands dominated 
by different friction mechanisms, i.e., when a new mechanism is "turned on" very quickly with 
the increase of frequency. Although highly hypothetical, such situations should not be written 
off completely. 

7 Discussion 

The shortness of Earth's day was undoubtedly beneficial for proliferation of biological life, 
making the daily temperature variation moderate. The situation may be drastically different 
for the growing class of potentially habitable super-Earth exoplanets. Due to the observational 
selection effect, the spectroscopically detected super-Earths are found mostly around lower- 
mass stars, whose habitable zones are inevitably narrower and closer. For such systems, any 
conclusion about potential habitability of a given exoplanet becomes especially uncertain, and 
the analysis becomes intricately involved with regard to such parameters as the amount of 
stellar irradiation, the chemical composition of a hypothetical atmosphere, and the internal 
heating. The rate of rotation is also a crucial parameter which for now remains unavailable 
from observation. The most advanced climatic simulations are based on a certain assumption 
of the spin-orbit state of the planet, e.g., a tidal synchronisation (1:1 resonance) is assumed, 
as in Selsis et al (2011). A tidally synchronised planet showing the same side to its host star 
has this side always exposed to plentiful irradiation, the other side always staying dark. Such 
planets can hardly retain an atmosphere and can hardly be habitable. However, a spin-orbit 
locking into higher commensurabilities (e.g., 3:2, as in the case of Mercury) allows the planet 
to rotate with respect to the host star, and leaves the planet a possibility of sustaining a stable 
atmosphere and having water in the liquid form on the surface. The 3D climatic simulations of 
the potentially habitable super-Earth GJ 581d, by Wordsworth et al. (2010), were performed 
for a set of possible spin-orbit resonances, including 2:1. This resonance was found to be the 
likeliest state of GJ 58 Id for a wide range of rheo logical parameters, assuming a terrestrial 
composition of its mantle (Makarov et al. 2012). Wordsworth et al. (2010) drew attention to 
the fact that a tidally synchronised atmosphere may be short-lived because of the collapse of 
CO2 on the night side. A similar phenomenon may occur on a slowly rotating planet. Both the 
constant angular lag tidal model and the constant time lag model predict that oblate planets 
with moderate and large eccentricities are captured in stable pseudosynchronous rotations, in 
which case their spin rate only slightly exceeds their orbital mean motion. This would make 
an entire class of detected exoplanets unsuitable for biological life. In this paper we prove that 
the prediction of pseudo-synchronism is germane to the above-mentioned simplified models of 
tidal interactions, models inapplicable to solid planets or moons. Super-Earth exoplanets of 
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Earth-like composition on eccentric orbits are likely to be captured into spin-orbit resonances 
higher than 1:1, but there is no such thing as pseudo-synchronous rotation for these objects. 
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Appendix A. The tidal torque 

The additional potential U of a tidally perturbed body can be expanded into a Fourier series 
over the tidal modes 

UJimpq = (I -2p) Co + (I -2p + q) M + m (SI - 9) « (I - 2p + q) n - m9 , (16) 

where 6 and 9 are the sidereal angle and rotation rate of the body, while oo , Q , n , and Ai 
are the periapse, the node, the mean motion, and the mean anomaly of the perturber as seen 
from the perturbed body. 

While the tidal modes oo lmpq can be of either sign, the forcing frequencies 

Xlmpq = | Ulmpq | ~ | (I ~ %P + <?) U ~ 771 9 | (17) 

at which the strain and stress oscillate are positive-definite. 

The series expansion of the additional potential U was developed by Kaula (1964), its 
partial sum known yet to Darwin (1879). Therefore the series for U and the resulting series 
for the torque are often named the Darwin- Kaula expansions. 

An accurate derivation of the expansion for the torque demonstrates that the torque con- 
tains both a rapidly oscillating and a secular part (Efroimsky 2012a). Having a zero orbital 
average, the oscillating part nevertheless may play a role in dissipation of free librations. In 
Makarov et al. (2012) it was explored whether the oscillating part of the torque can influ- 
ence capture into resonances. Changing the outcome of a particular realisation of the capture 
scenario, the oscillating part did not alter the statistics. 

The secular part of the polar torque looks as 

°° R l + 1 1 (I — m)\ 1 °° 
(%) = 2 G M* tar ~2T+2~ (i T m \\ m F imp( l ) Y G L( e ) H^impq) sin ei(u lmpq ) , (18) 

1=2 & m.=0 ^ p=0 q=-oo 

where the angular brackets signify orbital averaging, G denotes Newton's gravity constant, 
a, i, e are the semimajor axis, inclination (or obliquity), and eccentricity, while Fi mp {i) and 
Gi pq (e) are the inclination functions and eccentricity polynomials. The Love numbers k { and 
the phase lags e t depend on the modes ( |T6l) . 
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In the Darwin-Kaula theory, the phase lags come into being as products of the modes iOi mpq 
by the corresponding time delays: 

£l{Ul mp q) = (jJlm.pq ^Ui^lmpq) j (19) 

where, for causality reasons, the time lags At l (u l q ) are positive-definite. Therefore, (IT31) 
may be rewritten as 

€l{^lmpq) = Xlmpq A^li^lmpq) Sgn ( COl mp q 

) , (20) 

Ximpq being the physical forcing frequencies ffTTj) . As a result of this, the entire expression for 
the polar component of the torque can be written down as 

(%) = 



2 G^f s L Y + Y Y (i + m)'\ m ^ F *™p® Y G iv( e ) H u impq) sin | ei{u imp q) | Sgn ( u lmpq ) 

1=2 CL m=0 ^ p=0 q=-oo 

Be mindful that, similar to the Love numbers, we prefer to denote the lags with0 ei(uJi mpq ) 
and Ati(ui mpq ) , and not with ei mpq and Ati mpq . 

When the bodies are not too close (R/a <C 1), we drop the terms with I > 2 . For small 
obliquities (i ~ 0), we leave only %— independent terms (the next-order terms being quadratic 
in i). Finally, for eccentricities not exceeding ~ 0.3 , only the terms up to e 7 are important. 
Formally, this would imply summation over q = — 7, . . .7. However, the term with q = — 2 
vanishes identically, while those with q = — 7, ... — 3 are accompanied with extremely 
small numerical factors and can thus be dropped. So the polar component of the torque is 
approximated with 



3 7 

- GM s ] ar R 5 a-« G IJ & ) ^K 20g ) sin|e 2 (u; 220( p | Sgn ( ui 220q ) + 0(e 8 e) + 0(z 2 e) , 

<?=-! 

For the first time, this expression (with a sum running over all integer q = — oo, . . . oo ) was 
written, with no proof, by Goldreich & Peale (1966). A sketch of a proof was later suggested 
by Dobrovolskis (2007). 

The functional form of the dependence of the factors k\ sin e; upon the mode u: lrnvq is 
determined by the size and mass of the body and by its rheological properties. By rheology we 
imply the so-called constitutive equation of the medium, i.e., an equation interconnecting the 
strain and stress. For linear deformations, such equations can be rewritten in the frequency 
domain where each harmonic mode of the strain becomes expressed algebraically through the 

12 Although Kaula (1964) denoted the phase lags with e lmpq , the notation e/(w impg ) is more logical. It 
serves to emphasise the fact that for a homogeneous near-spherical body the functional dependence of a lag 
upon the tidal mode is defined by the degree I solely, while the dependence of the lag upon to, p, q comes 
about only due to the tidal mode w lmpq being dependent on these integers. 

While in the case of triaxial bodies the functional form of the lags is parameterised by all the four integers, 
for small triaxiality this complication may be ignored. 
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appropriate harmonic mode of the stress. Using the method explained in Efroimsky (2012a,b), 
the algebraic relations can be used to find the shape of the functions ki(u lmpq ) smei(u lmpq ) 
standing in the terms of the Darwin-Kaula expansion of the tidal torque. 

Calculations of the shapes of these functions, presented in Ibid., are based on a combined 
rheological model (Andrade at higher frequencies, Maxwell at lower frequencies), because of 
this model's ability to best match laboratory experiments and both seismic and geodetic mea- 
surements of dissipation over a broad range of frequencies in the solid Earth. It is reasonable to 
assume that the model is applicable to the mantles of other terrestrial planets and moons. As 
demonstrated in Ibid., this combined model furnishes for ki sine; a kink-shaped dependence 
upon the Fourier mode - as in Figure HI 



Appendix B. Calculation of the factors h(u lmpq ) s'mei(uJ lmpq ) 

As explained in Efroimsky (2012a,b), the products ki(u lmpq ) smei(u lmpq ) can be expressed 
through the mass and radius of the planet, and the real and imaginary parts of the complex 
compliance of its mantle. Thereby, the shape of the functional dependence of k\ sine; upon 
0J lmpq is defined by both the self-gravitation of the planet and its rheological properties. The 
functions turn out to be odd. They go continuously through nil and change their sign when 
the argument oo lmpq goes through nil, i.e., when a commensurability is crossed. 

These odd functions can then be written down as h(cu lmpq ) sin | €i(u lmpq ) | Sgn ( to lmpq ) ■ 
Here the product ki(u) lmpq ) sin | €[(u lmpq ) | is an even function of the tidal mode and can thus 
be treated as a function not of the tidal mode 0J lmpq but of its absolute value Xi mpq = I u im Pq I > 
which is the actual frequency of the oscillating stress in the mantle: 

h {u lmpq ) sin e z (u lmpq ) = h (u lmpq ) sin | e { (u lmpq ) | Sgn ( u lmpq ) 

= UXlmpq) s[n \ e l(Xlm P g) I S S n ( ^Im.pg ) ■ ( 23 ) 

The development in Ibid, results in the following frequency dependence: 



> . , v 3 -A, Jim [Jjx), 

where \ is a shortened notation for the frequency Xi mp q > while coefficients Ai are given by 



(2 I 2 + 41 + 3) 3 (2 I 2 + 41 + 3) /i 
IgpR " " 4 InGp 2 R 2 



A i = tz—f; - - A , _ Ti J2 d2 — • \ 2b > 



with G being the Newton gravitational constant, and R, p, p , g, and G being the radius, 
mean density, unrelaxed rigidity, and surface gravity of the planet. The functions 

- / cm \ 

Ke[J( X )} = J + J(xr A y a cos( — ) T(a + 1) (26) 

and 

lm[J( X )] = - J (xrj- 1 ~ J(xr A y a sin(y) T(a + 1) (27) 
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are the real and imaginary parts of the complex compliance J(x) of the mantle. Here a is 
the Andrade parameter assuming values of about 0.3 for solid silicates and about 0.14 — 0.2 
for partial melts. In our computations, we used a = 0.2 . The quantity J is the unrelaxed 
compliance of the mantle, which is the inverse of the mantle's unrelaxed rigidity fi . The pa- 
rameters t m and t a are typical timescale characterising the mantle's viscoelastic and inelastic 
response, correspondingly. 

The Maxwell time r u is the ratio of the mantle's viscosity r\ to its rigidity fi . While 
for the Earth's mantle it has a value of about 500 years, it may be much shorter for warmer 
planets and moons due to the exponential temperature-dependence of the viscosity. 

The inelastic (Andrade) time r A is expected to be of the same order as or lower than 
r M , over the frequencies higher than some threshold. For these frequencies then, the inelastic 
(containing r A ) terms in ( 1261 - |2T|) will be comparable to or larger than the viscoelastic (con- 
taining t m ) terms. However, at frequencies below the threshold, inelasticity ceases to play 
a major role in the internal friction, giving way to viscous friction which becomes dominant. 
Therefore at very low frequencies the mantle's behaviour approaches that of a Maxwell body. 
Mathematically, this implies that below the threshold the parameter t a increases rapidly as 
the frequency goes down (Efroimsky 2012a,b). So only the first term in ( 126]) and the first term 
in f )27|) are important, and we arrive at the complex compliance of a Maxwell material. The 
location of the frequency threshold may vary considerably for different planets. For the Earth, 
it is of the order of 1 yr _1 (Karato Sz Spetzler 1990). 

Numerical computations show that the probabilities of capture into resonances are not very 
sensitive to the value of r A , nor to the location of the threshold, nor to how quickly inelasticity 
yields to viscoelasticity with the decrease of frequency. In our numerics, we treat t a in the 
same way as in Makarov et al. (2012) and Makarov (2012). We set the threshold to be the 
same as for the solid Earth, 1 yr~ x . We then kept r A = t m over the frequencies above the 
threshold. For frequencies lower than the threshold, we set r A to grow exponentially with the 
decrease of the frequency. This way, at low frequencies the rheological model approaches the 
Maxwell one exponentially. For details, see Ibid. 

Writing a code, it is easier to divide both the numerator and denominator of (T2"4"j) by J 2 : 



3 —Ail 

Hu lmpq ) sine z (^ mp J = 2{l _ x) - ~ ~ ~ ~ 2 ~ j2 Sgn(w ImM ) , (28) 

where 1Z and X are the dimensionless real and imaginary parts of the complex compliance: 
K = 1 + ( X r A r a cos(^) r(a + 1) , (29) 

X = - (xrj- 1 - ( X r A )- a sin (^) r(a + 1) , (30) 
X being a short notation for the physical forcing frequency Xi mpq = I ^impq I ■ 
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